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We developed a gauge- covariant formulation of the non-equilibrium Green function 
method for the dynamical and/or non-uniform electromagnetic field by means of the de- 
formational quantization method. Such a formulation is realized by replacing the Moyal 
product in the so-called Wigner space by the star product, and facilitates the order-by- 
order calculation of a gauge-invariant observable in terms of the electromagnetic field. An 
application of this formalism to the linear response theory is discussed. 

§1. introduction 

The interaction between matter and electromagnetic fields constitutes the most 
important ingredients of condensed matter physics. Especially, since the high in- 
tensity fields become available, the non-linear and non-equilibrium responses are of 
great current interests. Therefore a microscopic quantum theory for these processes 
is called for. 

There already exist established theories for the response to external perturba- 
tions. The Kubo formalism 1 ) combined with the Matsubara Green function 2 ) is 
the most standard and successful linear response theory. However, its extension to 
higher-order responses to the non-equilibrium field is not straightforward. Another 
method depends on a semiclassical description such as the Boltzmann transport 
equation. 3 ^' 4 ) Although appealing to one's intuition, it sometimes leads to different 
results from what is obtained from the rigorous Kubo formalism. In contrast to 
these methods, the Keldysh formalism 5 )' 6 ) yields a powerful theoretical framework 
for directly handling the non-equilibrium Green function in the presence of electro- 
magnetic field. Kandanoff and Baym introduced the Wigner distribution function 
(WDF): f(X,p), where X^ = (T,X) and p^ = (u>,p) are center-of-mass coordi- 
nates and relative energy-momenta, respectively, and derived the equation of motion 
for the WDF. The WDF is an appropriate function to describe the interacting, 
many-particle system, and the equation for its time-evolution is called the Quantum 
Boltzmann Equation (QBE). However, they did not include the electromagnetic field 
in this formula. Altshuler, 7 ) Langreth, 8 ) Rammer, 9 ) and Mahan 10 )' 11 ) included the 
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electromagnetic field in such formula. However, it is not so straightforward to apply 
these formulas to the non- linear responses. 

Recently, we developed a generic and systematic theoretical framework for the 
linear and non-linear responses under the external constant electromagnetic field. 13 - 1 
In that work, we have employed the Wigner representations composed of the set 
of the center-of-mass time-space coordinates X and gauge covariant momentum 
7r = p — qA(X), where p is the momentum, q is the electric charge and A(X) is the 
electromagnetic potential. Then, we have accomplished the gauge-covariant formu- 
lation of the non-equilibrium Green function in the (X, n) space under the constant 
electromagnetic field, by replacing a conventional product with the Moyal product. 12 ) 
This framework gives a geometric view analogous to the general relativity: the ex- 
ternal electromagnetic field is incorporated in the geometry of the Wigner space. 
Corresponding to the commutation relationship [£^,7r^] = ih5u between operators, 
the Wigner space has the noncommutative geometry between time-space coordi- 
nate and energy-momentum. Furthermore, in the presence of the electromagnetic 
field, the non-zero commutator [itu, % v \ = ifoqF^ leads also to the noncommutativity 
within the energy-momentum space. It is highly desirable to generalize this geomet- 
rical formulation to the dynamical/non-uniform electromagnetic field. In this paper, 
we will show that this generalization is achieved by using a deformational quanti- 
zation method. 14 ^ In the linear order in the electromagnetic field, this formalism 
completely agrees with the Kubo formalism. ^ 

This paper is organized as follows. In £j2j we begin with a brief review of the 
Keldysh formulism and the gauge-covariant Dyson equations with external constant 
electromagnetic fields, which has been developed in the previous paper. 13 ) For self- 
containedness, a brief introduction to the star product and deformational quantiza- 
tion 14 ) is given in £j3j £0] constitutes the main body of the present paper, where we 
present the gauge-covariant formalism taking into account the generic electromag- 
netic field using the deformational quantization method. In an application to the 
linear response theory is given. 

§2. Keldysh formalism 

In this section, we briefly review the Wigner space of (X, ir) and the Moyal 
product in the case of the constant electromagnetic field, which have been developed 
in the previous paper. 13 ) 

First, we introduce a matrix of Green functions G in the Keldysh space, 

where the superscripts R, A, and < denote the retarded, the advanced and the lesser 
Green functions, respectively. They are defined as 

G R (x 1 ,x 2 ) ai , a2 ■.= -i6(t l -t 2 )([^ ai {x 1 )^i 2 {x 2 )] T ), (2-2) 

G A (x 1 ,x 2 ) ai ,a 2 ■=iO{t 1 -t 2 ){[ij ai {x 1 )^i 2 {x 2 )] T ), (2-3) 

G ,< (xi,x 2 ) Ql , a2 := Ti(i ) a 2 (x2)^>ai ( x i))i (2-4) 
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where ^ is a Bose (upper sign) or a Fermi (lower sign) field, ifv is a Hermitian conju- 
gate of tp, and x\ = (ti,xi), a\ represents an internal degree of freedom. Hereafter, 
we consider electron systems. The matrix G of the Green functions satisfies the 
Dyson equations: 



( £(0)-l 



£)*G 
I) 



ai,»2 



(x 1 ,x 2 ) = S aita2 5(x 1 - x 2 ), 
(xi,x 2 ) = S ai ,a 2 S(xi - x 2 ), 



(2-5) 
(2-6) 



a(°) 



where G is the unperturbed Green function for free electrons (note that the super- 
script (°) means the "unperturbed" in this paper), £_ is the matrix of the self-energies 
which arises from electron-electron and/or electron-phonon interactions and poten- 
tial scattering. (/ * g)(x\, x 2 ) := f dx^f{x\,x^)g{x^^x 2 ) is the convolution integral. 
Next, we introduce the Wigner representation: the center-of-mass coordinates 
= (T,X) := {\{t\ + t 2 ), |(cci + x 2 )), and the energy-momenta = (uj,p) := 
(uji —uj 2 ,p 1 —p 2 ) for relative coordinates. 9 )> 10 ) Then, the Dyson equations (|2-5p and 
(12-6|) are rewritten as 



m *hQ_) (x, P ) = i, 

(0)-l 



(G^(G W -Z))(X,p) = l, 
where the symbol represents the Moyal product defined by 

(/ *fi 9)ai, ct2 



(2-7) 
(2-8) 



(x, P ) = ^hu«Ax,p)e^ {9x ' 9p »- dp » 9x ' ) g a3 ,a 2 (x,p) 



"3 



(2-9) 



«3 



where d and d represent the derivatives which operate only to the left-hand side 
and the right-hand side, respectively. The average value of the physical observable 
which is represented by a one-particle operator can be calculated by a one-particle 
Green function. For example, the particle density p{X) and the charge current Jk(X) 
are written by G < (x) as 



p(X) = 
Jk(X) 



d d+1 p 
d d+l p 



tr 



G<(p,X) 



(2-kK) 



d+l 



tr 



V Pk H(p-qA(X)))G < (p,X) 



(240) 
(2-11) 



d x »A„ 



Under the uniform and static electromagnetic field F^ v = dx^A u — 
const., we have recently derived the gauge-covariant form of the Dyson equation 

(Gi 0) ~\x,Tr) -r(X,vr))exp(y( ( 9xM ( 9 7rM -t^x, + qF^^~d^))G(X,n) = 1, 

(2-12) 
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G(X,7r)exp(y(VxM?. M -t^~d x , + qF^t ^ \x, tt) - E(X, tt)) = 1, 

(2-13) 

where we have used the transformation of variables: 8 )' 13 ) 

Plt - qA^X) h-> 7r M , (2-14) 
f(X,p) ~ f(X,ir). (2-15) 

Here 7r^ is the mechanical momentum in the Wigner representation with the electric 
charge q of the particle, and the subscript o represents the zeroth-order terms with 
respect to F. All products appearing in the absence of the electromagnetic field 
is replaced by the Moyal product. 12 ) The operator exp(^( d d nfl — d nil d x>> + 
qF^ d ,rM d n *)) has been obtained from exp(y(<9 d Pfl — d p m d X n)) by chang- 
ing the coordinates from to ir^ and using fact that F^ v = const. However, with 
dynamical/non-uniform electromagnetic field, derivative d X i* generates dx^F^ 0) 
which complicates the expression of the product in the (X, tt) space. Below we will 
develop a rigorous treatment of these dynamical/non-uniform cases, which enables 
for example the description of non-linear responses to even oscillating electromag- 
netic fields. 

§3. Star product and Deformational quantization 

To achieve generalization to the case in general electromagnetic fields, we will 
use the deformational quantization method. 14 ) It provides the prescription to con- 
struct the star product from the usual product f-g and the Poisson bracket {/, g} := 
^ a' 3 ' (d x if) (d x j g) . Here, {/, <?} is assumed to satisfy Jacobi's identity: {/, {g, h}}+ 

ij 

{g, {h, /}} + {h, {/, g}} = and the matrix a* 3 is called the Poisson structure, 
where /, g and h are smooth functions on a finite-dimensional manifold M. Jacobi's 
identity can be rewritten by the condition for the Poisson structure: 

d 

(a a d x ia ik + a jl d x ia kl + a kl d x ia ij ^j = 0, for k = 1 ~ d, (3-1) 

l=i 

where we have used a* 3 = —a 3 *, and d is the dimension of the manifold. The star 
product is written as 

oo 

f*g:= Y,C n (f,g)(ih/2r, (3-2) 

n=0 

where C Q (f,g) := f ■ g and C\(f,g) := {/, g}. Then C n (f,g) with n > 2 is defined 
so as to satisfy the associatively property: 



f*(g*h) = (f*g)*h. 



(3-3) 
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For our purpose, the manifold M is the Wigner space. Kontsevich 14 ) found an 
ingenious method to construct the star product, which we introduce here. 



Kontsevich's star product 

Kontsevich invented the diagram technique to calculate the n-th order term in 
ih/2 for the expansion of / ★ g: 



(f*g)(x)=f(x)g(x) + J2 (t) E vrBrM9), 



n=l 



ih 



(34) 



r&Gn 



where r, Br t a(f,g) and wp are defined as follows: 14 ) 
Definition. 1 

G n is a set of the graphs r which have n + 2 vertices and 2n edges. Vertices are 
labeled by symbols "l",---,"n", "L", "R". Edges are labeled by symbols (k,v), 
where k = 1, ■ • • ,n , v = 1, ■ • ■ ,n,L,R and k ^ v. (k,v) means the edge which 
starts at "/c" and ends at 'V. There are two edges starting from each vertex with 
k = 1, ■ • • ,n. L and R are the exception, i.e., they act only as the end points of 
the edges. Hereafter, Vp and Ep represent the set of the vertices and the edges, 
respectively. 

Definition. 2 

Br,a{f->9) IS th e operator defined as: 



B r>a (f,9):= 

I:E r ^{i u - 



n n ■ 

k=l \e£E r ,e=(k,*) 



71(e) 



a 



I((k,vi),(k,vD) 



n 



7/(e) 



/ 



yee£r,e=(*,L) 



II d m ) a 

K e&E r ,e=(*,B) 



(3-5) 



where, / is a map from the list of edges ((k, v^ 2 )), k = 1, 2, • • • n to integer numbers 
*2j • * * 5 *2n}- Here, 1 < i n < d and d is a dimension of the manifold M. Br, a (f,g) 
corresponds to the graph r in the following way. The vertices "1",- • • , "n", corre- 
spond to the Poisson structure a ij . R and L correspond to the functions / and g, 
respectively. The edge e = (k, v) represents the derivative acting on the vertex v. 
Thus the graph r corresponds to the bidifferetial operator Br, a (f,g) by the above 
rules. The simplest diagram for n = 1 is shown in Fig. [1] (a), which corresponds to 
the Poisson bracket {/, g} = a %xl2 {d x % 1 f)(d x % 2 g). The higher order terms are the 

generalizations of this Poisson bracket satisfying Eq. (|3 3|) . Fig. [2(b) shown a graph 
r with n = 2 with the list of edges 



((1,L), (1,12), (2,12), (2,3)). 
The operator Bp t a corresponding to this graph is 



(3-6) 
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(a) 





Fig. 1. (a): The graph fed corresponding to Poisson bracket, (b): A graph r £ G2 correspond 
to the list of edges: ((1, L), (1, R), (2, R), (2, 1)) ^ {ii,ia,is,U}- 

(/,«/) ~ E (^^^(d^f^d^d^g). (3-7) 

Definition. 3 

We put the vertices in the upper-half complex plane H + := {z E C|Im(z) > 0}. i? 
and L are put at and 1, respectively. We associate a weight wp with each graph 
r E G n as 



1 



where <ft is defined by 



Lh 1 T __ f (q-p)(q-p) 



p and g are the coordinates of the vertices "fc" and 'V , respectively, p is complex 
conjugate of p. 7^ n represents the space of configurations of n numbered pair- wise 
distinct points on H + : 

ftn ■= {(pi, ■ ■ ■ ,Pn)\Pk G H + , p k / pi for h ± I). (3-10) 

Here we assume that H + has the Poincare metric: ds 2 = (d(Re(p)) 2 +d(lm(p)) 2 )/(lm(p)) 2 , 
p G i2" + . From the metric, we can calculate a geodesic. Note that, ^(p, q) is the angle 
with the geodesies which is defined by (p,q) and (00, p). Namely, <f) h (p, q) = Zpqoc. 
For example, wr corresponding to the Fig. [1] (a) is 



"> = w 1 4 L °« (?) A 4 L °« ( 1 = (3 - u) 



where we have included the factor "2" arisen from the interchange between two edges 
in r. wr corresponding to the Fig. [1] (b) is 

1 f f pi\ a j 1 t A-r^ , A /M . ,1, ( {pi - p2){pi - P2) 

w r (b) = n „ n u / d-Log — A d-Log — A cZ-Log — A d— Log 1 



(b) 2!(2vr) 4 J H2 i \p 1 ) i \l-p 1 J i \p 2 J 2i \(pi - P 2 )(Pi ~ P2) 
a-Log — A a-Lor 



2!(2vr)4y w2 i a \pJ i °V1-Pi 



Gauge covariant formulation of Wigner representation through deformational quantization! 
Ad(2arg(p 2 )) A d\p 2 \ — — r7TT L °g 



d\p 2 \2i \(p 1 -p 2 )(Pi ~P 2 ] 
1 d-Log ( ^ J A d-Log ( — — ^ ) A ci-Log ( ^ ) A d-Log ' ' ~ P ' 2 



2!(2vr) 4 J H2 i \p x ) i \l-p 1 J i \p 2 J i \l~P 2 

= "21 
1 

~ 2' 

where pi and £>2 are the coordinates of vertices "1" and "2" , respectively. Here, we 
have used the following facts: 

(Pi -P2)(Pl / (pi - ^^fe))^ _ yl e iarg( P2 )) 



(3-12) 



P. ,„ t f (Pi ~ P2)(P! ~ P2) \ j 

/ d p2 0|p2|Log 7= -p — r = hm Log 

Jo \\Pi ~ P2AP1 ~ P 2 ) J A-00 



(pi - yle-^s^ 2 ))^ - yle- iar s(P 2 )) 
(1 — 7le iarg ( p2 ))(l — yle iarg ( p2 )) \ 



— 1.1 j 11 IjOS! 

yl-+oo I (1 - yle- iar s(P 2 ))(l - yle-* ar s(P 2 )) J ' 

dLog f ^ A dLog ( ^ [ dLog ( A dLog ( = 0. (3-13) 
| P1 |>/1 VPi/ \1-Pi/ ^bi|>yl VPl/ 

Generally speaking, the integrals are entangled for n > 3 graphs, and the weight 
of these are not so easy to evaluate as Eq. (|3-12|) . 

§4. Gauge-covariant formulation with electromagnetic fields 

In this section, we derive the gauge-covariant Dyson equation. Since we can use 
the deformational quantization method, we only need Cq and C\ in (|3 2j) . We should 
calculate the zeroth-order and the first-order terms which are given by expanding 
the Dyson equation (|2 6j) with respect to h and by the variable transformation p^ — 
qA^(X) 1— > 7r^. The derivation is organized by three steps as shown in FigO 

4.1. Transformation of variables and Poisson structure 



Step. 1: We derive the zeroth-order and the first-order terms of the Dyson equations 
with respect to h by the gradient expansion. Expanding the left hand side (LHS) of 
Eq. (|2-5p with respect to h, we obtain the first two terms as 

((G^~\x)-l(x))-G(x)), (4-1) 

|{(g!, 0K1 W-^(x)),G(x))}, (4-2) 

where x := (X,p) and the dimension of the Wigner space is 4 + 4. {/, g} := 
^2 a o (9 x if)(d x jg) is the Poisson bracket with the Poisson structure 0$ = 

1 ~ 4 + 4) given by: 
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(X, P ) 



(X,*) 



f{X,p)-g{X,p) 

{f,g}(x, P ) 



^> !/-;y|: v..-; 



/(A»- 9 (X,7r) 



(f*ng)(X, P ) 



(J*g)(.X,n) 



Fig. 2. This figure represents the derivation of the LHS of the gauge-covariant Dyson equation, 
where / and g are certain functions in the Dyson equation. The derivation has three steps. 
First, we obtain the Poisson structure of the Dyson equation with Wigner space (X,p). This 
Poisson structure is obtained by using the gradient expansion, 9 ' which is represented by the 
symbol "-f|-". Secondly, we perform the variable transformation, p M — qA^(X) i— » 7r M , which is 
represented by the symbol . Thirdly, explicitly gauge-covariant Dyson equations are given 

by using the deformational quantization method which is represented by the symbol "JJ." . 

and if is the Minkowski metric i) = diag{— 1, 1, 1, 1}. 

Step. 2: Next, we introduce the transformation of variable f)2 • 14|) and (|2- 15[) . This 
step corresponds to the calculation of Cq and C\ in definition f|3-2|) in (X, tt) space. 
Using the transformations (|2-14|) and (|2-15|) . we can write the derivative of / and g 
as follows: 



(d x ,f(X,p))(d^g(X,p)) - (d p J(X,p))(d x ,g(X,p)) 
= (d x »f(X,7r))(d^g(X,7r)) - (d n J(X,7T))(dx,g(X,7r)) 

-q ((d x »A u (X))(d„J(X,iT))(d^g(X,TT)) - (^^(X))(a^/(X,7r))(^ 5 (X,vr))) 
= (dx,f(X,7T))(d^g(X,7T)) - {dvJ{X^))(dx»g(X^))+qF^J{X,i:)){d^g{X^)), 



where F is the field strength: = dx^A u — dx^A^. Accordingly the Poisson 
structure is transformed as: 



By this replacements, the zeroth-order terms and the first-order terms are trans- 
formed as 



(4-4) 





qFf 




(4-5) 



*(0)-l 



(5)-£(x))-G(x)) 



(4-6) 
(4-7) 
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where x := (X,tt). In (|4-7j) . {/, g} := a 1 ^ {d 5 , 1 f){d 5 , ] g) is also the Poisson bracket, 

i.e., satisfies Jacobi's identity. Hereafter, the tilde in x will be omitted for simplicity. 
4.2. Gauge- covariant formulation 

Step 3: From the Oth-order and the first-order terms obtained in the previous section, 
gauge-covariant-Dyson equations are given by using the deformational quantization 
method. By constructing the star product from (|4-6p and (|4- 7|) . the Dyson equation 
reads: 

m^~\x,n) -Z(X,-k))*G(X,tt) = 1, (4-8) 

where the symbol is star product which is given the Poisson structure a with 
the Wigner space (X, it). 
Similarly, 

G(X,ir)*((G^~\x,n) -£(X,n)) = 1. (4-9) 

Note that the difference between Eqs. (|2-7p and (|4-8p is only in the definition of 
the product. Namely, to obtain the explicitly gauge-covariant Dyson equation, we 
only change the product to 

Response to electromagnetic field 

Let's recall the diagram rules for Kontsevich's star product. The rules relates 
a bidifferential operator to a graph r € G n . G n is a set of nth-order graphs with 
respect to h, and each graph r G G n is constructed by n+2 vertices which are labeled 
as "1", "2", "n", "L" and "R", and by 2n edges which start from the vertex 
u k" , k L or R. The vertices which are labeled as "1" , "2" , • • • , "n" , correspond to 

a. The vertices which are labeled as "L" and "i?" correspond to G or (Gg°^ — IT). 
We separate the Poisson structure into two parts ao + of as (|4-5j) and expand G 

and (Gg°^ — in terms of F: 

G = Gq + G p + G pi + • • • 

-Z) = - A)) - If - ^ ~ ■ ■ ■ , (4-10) 

where subscript o and F k represent zeroth order and A:-th order terms with respect 
to F, respectively. Thus the contributions corresponding to each graph are further 
classified according to the order in F. Now, the vertices "1" , "2" , ■ • • , "n" correspond 
to ao or aF- The numbers of and «o ar e denoted by n ap and n ao , the vertices 
"L" and "i?" correspond to f F i and 5F m ; receptivity, where / and g represent G or 

G=lq°^ — ^), and l,m are some nonnegative integers. The order of the graph with 
respect to F is given by I + m + n aF . The Kontsevich's diagram rules applied to the 
present case read as 

1. n = n ap + n ao vertices "1", "2", • • • , "n" are put in the upper-half complex 
plan H + := {z € C|Im(£) > 0} and two vertices "L" and "i?" are put at and 
1, respectively. 
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2. Vertices "1", "2", V have two edges (k,vl' 2 ),k = 1 ~ n,vl' 2 = l~ n ,L,R 
and k ^ v\,k ^ v"l, where each edge (k, v) starts at vertex "fe" and ends at any 
other vertices V . 

3. The set of vertices and edges are written by r. G n is defined by the set of all 
graph of n + 2 vertices and 2n edges. G n G r. 

4. Vertices "1", "2", • • • , "n QF " correspond to a F , vertices u n ap + 1", "n QF + 2", 
■ ■ ■ > " n a F + n "o = n " correspond to ao, and "L" and "i?" correspond to f F i 
and <7_F m ; respectively. 

5. The operator B Fa is constructed as follows. Edge (k, v) corresponds to the 
differentiation with respect to x acting on a function corresponding to vertex 
V (I = 1 ~ n,l ^ k), "L" or "R". Thus a set of the vertex "A;" and the two 

edges (/c,f^' 2 ) corresponds to bidifferential operator ^ a Q OT p9 x i(h v i )d x3 >(h v 2), 

ij k k 

where h v is the function corresponding to the vertex v. 

6. The weight wp of r is defined by 

wp '■ 



" ^ Log V(5-p)( 5 -p)J ' 



where p and q are coordinates of vertices "A;" and "u" in H+, respectively, and 
H n is the space of configurations on n numbered pair-wise distinct points on 
H+: H n = {(pi,-- - ,Pn)\Pk G H + ,p k ^pi for k ^ /}. 
7. The rth-order star product in terms of F is given by 



oo 



(f Fl *g Fm ) F r= (. ih / 2 ) naF+na ° Yl W r B F:a (f Fh g F m), 



where r = I + m + n QF . 



Additional diagram rules: 

Al. Two edges of each vertex corresponding to connect with both vertices "L" 
and "12". 

A2. Edges of the vertices corresponding to qo do not connect with the vertices 

corresponding to ao, because r] is the constant metric. 
A3. At least one edge of the vertex corresponding to ao connects with "L" or "R" . 

To calculate (/ * g) F r \ /=/ , , first we separate the graph r into r ap and r ao . 

g—gpm 

r ap is the graph consisted by vertices corresponding to aj?, and "L" and U R" , and 
edges starting from these vertices. r ao is the rest of the graph r without r ap . 
Namely, r a „ is the set of the vertices corresponding to «o which are labeled by "fe" 

12 

(k = n ap ~ (n ap + n Qo )) and edges (k,v^ ). Next, we calculate the weight w nap 
and the operator B Fa ^ ap corresponding to r ap , and later those for r ao . 



Separation of graph r 
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2 
-A 



9F- 



Fig. 3. A four vertices graph, where the white circle and the white triangle represent ao and chf, 
respectively, and the dotted arrow and the real arrow represent the derivative with respect to 7r 
and X, respectively. 



We now prove w r B r , a = w naot n aF Br ao , ao ■ w naF Br aF , aF , where w„ Q0 ,n QF = 

J IT-CtQ 71 <^ J? 

(nJ+n Qo) ! • 5^T» w n F = and Wl is given by Eq. d^H]). 

First we consider the case where the edges of the vertices corresponding to ao do 

n \n ! ™ a ° 

not connect with a f- In this case, w n „ n „ = -, — °°' 1 , . Secondly we consider 

' n a ,n aF (n aF +n ao y. n aQ \ ■> 

the graph which consists of four vertices corresponding to ao, ctF, and "L"and "i?" 
as shown in Fig. [3J We also assume that one edge of the vertex corresponding to 
ao connect with vertex corresponding to ap. In this case, from additional diagram 
rule A3, another edge of the vertex has to connect with "L" or "i?". Since we can 
exchange the role "i?" and "L" by the variable transformation p *— > 1 — p, 
we can assume that one edge of the vertex corresponding to ao connect with "L". 
The weight in this case is given by Eq. (|3-12p . i.e., the integrals for the weight 
is given by replacing coordinate of the vertex corresponding to atF by coordinate 
of "it!" (or "L") in H + . From additional diagram rules, above result of integrals is 
applied to every graph. Namely, the position of each vertex corresponding to ao 
and cxf can be move independently in integrals, and the entangled integral does not 
appear. Therefore the weight w ntXQjna of a graph r ao only depend on the number 
of vertices corresponding to ao and of- and wr = Wn aQ ,n aF ■ w naF holds generally. 

Finally, we can obtain w Unn = > ap ' J , and separate the operator wpBp a 
into w nao B raoiao and w naF Br aF , aF . 

Explicit expression for star product 

We will proceed to obtain the main result of this paper (Eq. (|4-16[) ) for the gauge- 
covariant star product. We treat of, fpi and gp m as a cluster, and ao as an operator 
acting on the cluster as shown in Fig. HI In this figure, the operator corresponding 
to the graph (a) is calculated by two steps. First, the graph (6) corresponding to 
w n aF Br olF ,a F is treated as a cluster, and is given by Kontsevich's rules. Secondly, 
we calculate the operator Wn aQ Br ao ,a that acts on the cluster as in graph (c), where 
the big open circle represents the graph (b). Namely, 

(ih\ n 
-w ) w r ao Br ao , ao ■ wr aF Br ap ,«, (f F i ,9f^), 

"-an an ^ ^ T 
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(a) ») (c) 



2 




Fig. 4. This figure shows the calculation method of the graph (a), where the dotted arrow and real 
arrow represent the derivative with respect to n and X, respectively, and the white circle and 
the white triangle represent Qo and of, respectively. We rewrite the graph (a) as the graph (c) 
which is given by the cluster represented by the big circle and the operators into it, where the 
big circle represents the graph (6). 



(4-11) 

where G' is the set of r ao . 

We consider the contribution of the cluster. Since edges of ap(X) correspond 
to the derivative with respect to 7r^, Qi?-vertices do not connect with each other. In 
this case, the operator w naF Br apj a F is given by 



2 \ na F i ih<— — > 



(4-12) 



Using Eq. (|4-12|), (f*g) F r \ f =f pl is given by 

9— gpm 



— ) w naotnaF Br ao<ao ■ w naF Br aF ,a F {fFh9F" 



y-^ (n ao + n aF )\ / ih\ a ° " F n ao \n aF \ a Q c ^ ^ ^ ^ ^ 
Tl aQ \n aF \ y 2 J (j^cta ^af )■ ^ao' ' 



exp( ^a ) {^hi^^qF^^T^gpr, 



(4-13) 



where chq is defined by 



ao(/i, h,---)-= ^2 ^2 a oh--- {9 x ifki) ■ ■ ■ {d x jfk 2 ] 

l<fcl<fc2<oo i,j 



(4-14) 



Thus the rth-order terms (/ * g) ft is given by 

{f*g)Fr= Yl (f * 9)ft\ f=f Ft ■ (4-15) 



;,mSN 
0<i + m<r 



g — gpm 



Gauge covariant formulation of Wigner representation through deformational quantization 



Substituting Eqs. (I4-I3j) and (I4-I5j) into Eqs. (H^l) and fljjQ]) , we obtain the gauge- 
covariant Dyson equations in the presence of F^ v as 

(/ * g)(X, it) = e% & ° (f(X, ir)ef yr)). (4-16) 

ih * 

When /, g and i* 1 are expanded by Fourier series, the operator e^~ a ° equals to the 
translation operator in the momentum space. We use this fact in the following 
section. 

§5. Linear response theory 

As an application of the formalism developed above, we discuss in this section 
the linear response to the electromagnetic field F^ v oc e 1 *^^ starting from the 
Dyson equation. The linear components of Eq. (|4-8p read 

((Cj 0-1 -Z)*G) F (X,n)=0, (5-1) 
(G*^" 1 -Z))f(X,k) = 0, (5-2) 

where (tt) = it — H(tt) , H is the Hamiltonian of electrons and the linear order 

terms in can be written as 

(/ * 9)f ■= Uf *o go) + (/o *o 9f) + (fo *f go), (5-3) 
Each term is defined by 

(f F *ogo) := ef &0 (f F (X,7T)g (X,7T)), (5-4) 
(/o*o2f) :=ef &0 (f (X,TT)g F (X,7T)l (5-5) 

(fo* F go) := ef *°(f (X,ir)(j *d^qF^^)g {X,ir)), (5-6) 

where / and g represent G or (Gq^ Hereafter, we assume that the equilibrium 

system is uniform, and the nonequilibrium one is steady, i.e., f F (X, it), g F (X, tt) oc 
e iQ M x^ / (X,tt) = / (tt) and g (X,ir) = g (ir). In this case, Eqs. ([5% and 
(|5-6|) turn into 

(fF*ogo) = fF(X,TT)g (ir-Q/2), (5-7) 
(/o*o5f) = f (TT + Q/2)g F (X,7T), (5-8) 

(/o *F So) = y/o(vr + Q/2)*d ^qF^~d „ u g Q (TT - Q/2), (5-9) 

where e a ^^ in e~2~°° shifts the momentum 7r^ to tt^ + a^. From Eq. (|5-ip . the 
equation of the Green functions is given by 

G F (X, tt) = G^tt + Q/2)t F (X, tt)G (tt - Q/2) 

-G (tt + Q/2)((G (tt + Q/2)*d^qF^~d nv G (7T - Q/2)). 

(5-10) 
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We note that Eq. (15-lOp is rewritten by 




where riFermi is the Fermi distribution function and Gq = (G A — G^)riFermi- We have 
used the retarded and the advanced components of Eqs. (15-ip and fl5- 2j) : 



(Go^^)- 1 *o Gp ,A + (Gp' A )~ l * G^ A + (Gq' A )~ 1 * f G*> A = 0, (5-12) 
Gp' A * (G^^)- 1 + G^ A * {gVY 1 + G^ A * F (d**)- 1 = 0. (5-13) 



Eq. (|5-18p is the generalization of the Stfeda formula. 13 - 1 ' 15 ) 
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